Probability distributions in statistical ensembles with conserved charges 
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Macroscopic data of an equilibrium state are described by means of statistical distributions of microscopic 
variables specific for a given ensemble. In the application to the description of particle production in high energy 
particle collisions we are generally dealing with the grand canonical ensemble with respect to the particle number 
HI H Q ■ In the ultrarelativistic situation energy conservation and particle number are usually controlled by the 
temperature of the system Q . 

Applications of statistical physics concepts to multiparticle production processes require however the imple- 
mentation of internal symmetries [5|,|6| that result in associated conservation of quantum numbers. In the grand 
canonical formulation (GC) conservation of quantum numbers is implemented on the average and is determined 
by the corresponding chemical potential. On the other hand in the canonical ensemble (C) quantum number 
r"| , conservation is treated exactly. 

In general the thermodynamic quantities calculated in the GC and C ensembles differ. This is particularly the 
case when dealing with small systems |l|, |2|, y, |4J, []J |8| • Since relativistic heavy ion collisions correspond to a finite 



volume and to a given charge value, the canonical ensemble should be used whenever the GC and C formalisms 



give different answers [1 H M fiol |llj . It is however to be expected that in the thermodynamic limit the C and 
GC descriptions should provide the same answer for physical observables |2l vi llfjj. The thermodynamic limit 
is reached in the large volume for fixed density limit in C and for fixed average density in the GC ensemble. 
Only in this limit can one indiscriminately use the GC or C descriptions. The results presented here do not 
rely on the system under consideration being of a relativistic nature and apply equally well to relativistic and 
non-relativistic systems. 

There is an essential difference in the volume dependence of observables in the GC and C formulations 
HIl H 110 Consequently in the limit when V — > oo some ratios of extensive quantities could in general 
converge to different values in GC and C ensemble • It is thus clear that the equivalence of both descriptions 
in the thermodynamic limit can be strictly established only for intensive observables 0. 

In applications of statistical models to particle production in high energy collisions of elementary particles 
|13| and in heavy ion collisions P, IT3 | we are always dealing with small systems. Thus, the model description 
of particle yields is in principle canonical with respect to the conservation laws. This is particularly evident in 
elementary particle and peripheral heavy ion collisions where e.g. strangeness production is strongly suppressed 
due to canonical effects ^l). However, a detailed analysis of different particle yields 0,H EH EH El nas shown 
that in central heavy ion collisions the relative error between the C and GC descriptions is so small that the 
GC approximation can be used with confidence. 
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Recently, an interesting observation has been made [12| that in finite systems even very small relative errors 
between C and GC results seen on the level of first particle moments i.e. thermal particle multiplicities do not 
necessarily guarantee that this is also the case for higher moments or particle fluctuations. Contrary, for a set of 
thermal parameters relevant in high energy heavy ion collisions the second moment differs substantially in the 
GC and C ensembles at finite V even though the first moments agree to high accuracy. The above observation 
has been generalized ^2 to the large volume limit with the conclusion of the violation of GC and C ensemble 
equivalence in the thermodynamic limit on the level of particle number fluctuations. 

However, a direct comparison of GC and C results for the average charged particle number or higher moments 
is only adequate for a finite system. In the thermodynamic limit such a comparison can be only done at the 
level of densities. This is simply because in the large volume limit only particle densities are finite in GC and 
C ensemble whereas corresponding multiplicities are infinite in both cases. In the thermodynamic limit the 
equivalence problem of GC and C formulation of the conservation laws requires intensive observables. 

We show that different statistical ensembles with respect to conservation laws are exactly equivalent in the 
thermodynamic limit. Consequently, there is equivalence of all possible moments, relative fluctuations, scaled 
variance etc. The exact calculations are performed in the GC and C ensemble with arbitrary values of the 
conserved charge. 



II. PROBABILITY DISTRIBUTIONS IN CANONICAL AND GRAND CANONICAL ENSEMBLES 



The equivalence of canonical and grand canonical ensembles of statistical systems with conservation laws 
will be discussed for a non-interacting relativistic gas of charged particles and antiparticles in a volume V at 
temperature T. Particles and antiparticles have charge ±1 respectively. 

Let us consider a system with total overall charge Q with N = N- antiparticles and N + = N + Q particles. 
The requirement of the exact charge conservation in the canonical ensemble leads to the following partition 
function 0,0,13 

00 7 2N+Q 

^T)=Tr Q e-^= £ ^=1,(2,), (1) 
where z is the one-particle partition function 

AT) = -£pJ*p °~ 0V ^ - h Tm ^ (?) - v*oF) ■ ( 2 ) 

Iq and K2 are the modified Bessel functions 0] . 

To obtain the average particle multiplicity moments one introduces in e.g.Q) fugacity parameters A + and A_, 
which are set to one in the final formulae 

00 y 2N+Q / \ \ Q/2 

z-(v,t,x.,x +) = £ ^rNTm xN + +QxN = (r) *><WTO, (3) 



such that 

1 / a \ k _ 

(4) 



WI)q = ^b (a±^) -Z§(V,T,A_,A + ) 



From Eqs. (00} it is clear that 

J2N+Q i 

^N,N + Q,V) = m{N + Q)UQ{2zy (5) 

is the probability distribution to find N negatively and N + Q positively charged particles in a system of volume 
V, temperature T and total charge Q Hi. 17 . 

In the grand canonical ensemble the charge is conserved on the average, thus the partition function 

Z GC (V,T) = Tre- p6+ T3, (6) 
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where 7 is chosen such as to reproduce the average charge (Q) in the system. For a non-interacting gas the 
trace in the above equation can be calculated explicitly yielding 

00 00 „7(AT + -iV_) \ N + \ N - 

Z GC {V,T n ,X-,K) = E E ~ N w , Z ^ +N+ = exp[(A + e-+A_e^)z] (7) 

AT + =OAT_=0 + ' 

where as in Eq. ((3J) we have introduced auxiliary fugacities A± for particles and antiparticlcs. Following Eq. 
the average number of particles and the average charge in the GC ensemble is obtained as 



(N±)=z exp (±7), (Q) = (N+)- (iV_) = 2z cosh 7 
In terms of the total charge Q = N + — N- the GC partition function J7J) can be written as 



00 00 



Z^(T/T,7,A_,A + )=E E 



Nl(Q + N)\ 



N=0 Q=-N 

with iV = iV_ being the number of negatively charged particles. Thus, the function 

1 



^,N + Q,V)- ZGC 



7 Q+2N 

e~< Q 



(8) 



(9) 



(10) 



defines the probability distribution in a GC ensemble, with average charge (Q), to find a charge Q with N 
negatively charged particles. Expressing the chemical potential appearing in (|10Jl through the corresponding 
average charge from JSJ one finds 



r 2N+Q 



V$(N,N + Q,V) 



N\(N + Q)\ 



1 Q 



2z 



(ii) 



All other less restricted probabilities can now be obtained directly from Eq. 1)11(1 . The probability distribution 
(N, V) to find N particles or V^(Q, V) to find charge Q in the volume V at a given average charge (Q) 
are obtained from 

00 

V^(N,V) = E Vf^{N,N + Q,V), 



-PfQ){Q,v) 



E V%f } (N,N + Q,V). 

7V=max(-Q, 0) 



(12a) 
(12b) 



The summations in Eqs. (|12a|) and l|12bj) can be done explicitly yielding 

-1 N± 



For the charge distribution one finds 



V(Q) 2 +4z 2 ±(Q) 



exp 



v/(Q) 2 + 4z 2 ±(Q) 



Vffi (Q,V) = I Q (2z) 



iQ 



2z 



V<Q> 2 + 4z2 



(13a) 



(13b) 



The particle number probability distribution l|13al) is, as expected [Til Il7j. a Poisson distribution. The charge 
distribution, on the other hand, is not Poissonian due to the constraints imposed by the requirement of the 
exact charge conservation in a given sector of GC ensemble with a fixed average charge (Q). 



III. PROBABILITY DISTRIBUTIONS IN THE THERMODYNAMIC LIMIT 



In the previous section we have introduced a set of probability distributions in the GC and C ensembles 
related to the charged particle number. All these distributions are also valid for any value of the conserved 
charge as well as the volume of the system. Thus, they can be used to analyze the thermodynamic limit. 
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From (J5J and one finds 

Pffi (N, N + Q,V)= V$(N, N + Q,V) Vf^ (Q, V) . (14) 

It is straightforward that in the sector of fixed charge Q the corresponding particle number is distributed as in 
the canonical ensemble. That is why in Eq. (|14fl the GC probability function is just the product of the 
canonical particle number and the grand canonical charge distribution. 

To take the thermodynamic limit in l|14|) one first expresses the variables (N,Q, (Q)) by means of the cor- 
responding densities (n,q, (q)) and then one takes the limit V — > oo for fixed densities. This also requires the 
replacement of a discrete sum (1/V) J2n ~* J 

The essential difference between GC and C distributions in Eq. (|14fl appears through the probability function 
P/q? (Q, V). Thus, to study the equivalence of GC and C ensemble on the probability level in the thermodynamic 

limit, it is sufficient to find VfjF (q) from 

Vf q f{q) = Jiui VVf q f{q,V), (15) 

where is obtained from Eq. (|13b|) by replacing (Q) and Q by V(q) and Vq respectively. The large volume 
limit in 115|) is taken at fixed densities. An extra volume factor in Eq. (|15|) appears from the replacement of 
discrete by continuum variables. 

The limit in Eq. (|15fl is obtained from the a — > oo behavior of the Bessel function pj| 



I a (ax) 



gCt\/l+a; 2 



Consequently, the charge density probability distribution 
V? Q ? (?)= lim IV 1 ' 2 * 



i + vT+~ 



V-»oo | % /2T(q 2 + 4^) 1 /4 



q+^q 2 + 4z 2 



Vq' 



(16) 



(17) 



It is rather straightforward to see from Eq. I|17|) that the limit V — * oo does not exist as a regular function. 
This limit is zero for any q ^ (q) and infinity for q = (q). Let us consider however the thermodynamic limit in 
Eqs. (|15|l and l|17(l as a generalized function 

T = K fim V J dqpVf{q, V)f(q) . (18) 

The density integration in l|18|) is obtained through the saddle-point method 



T = Mm V 1/2 / dq e vs(q) /(g), (19) 



where the function 



S(q) = ^q 2 +Az 2 -^(q) 2 + Az 2 + q log ({q) + tJ (q) 2 + 4z 2 ^ - q log ^ + ^ q 2 + Az 2 }j 
In the large volume limit the dominant contribution to the integral l|19|) is obtained as 



(20) 



* = Jte. 1,1/2 { f vfe 3{ -"> eVS< " > +0 ^\- < 21 > 

where q = (q) is just a saddle-point such that S (go) = 0. From Eq. (|^T|) one finds 

^ =/««)), 

which means that the charge density probability distribution (|15|l converges to a delta function 

rf q r(q) = S(q-{q)), (22) 

The above result together with Eq. (|14|) taken in the thermodynamic limit completes the proof of equivalence 
of the GC and C ensembles on the probability level. The probability to find a given density of particles and 
antiparticles in the GC ensemble with a fixed average charge density (q) is exactly equal to the corresponding 
probability in the C ensemble if one identifies charge density q in the C ensemble with (q) . The same replacement 
q — > (q) is required for all intensive thermodynamic observables. 
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IV. SUMMARY 

The equivalence of the grand canonical and canonical descriptions of the conservation laws has been considered 
in the thermodynamic limit. The problem has been studied in the relativistic gas composed only from one type 
of non-interacting particles and antiparticles with the constraints imposed by charge conservation. Detailed 
studies of the equivalence problem of canonical and grand canonical ensemble were presented on the level of 
probability distributions. In the thermodynamic limit the probability distributions for particle and antiparticle 
densities coincide in canonical and grand canonical ensembles. Consequently, for charged particle densities there 
is full equivalence of all possible scaled moments, relative fluctuations and scaled variance in the thermodynamic 
limit. 
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